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Abstract. We report the detection of ~1 500 RR Lyrae 
of Bailey type ab located in the Sagittarius dwarf galaxy 
(Sgr). These variables have been detected on two ESO 
Schmidt fields centred on (l,b)=(3.1° r 7.1°) and (6.6°,- 
10.8°), covering an area of ~50 deg 2 . We present a surface 
density map of Sgr based on the spatial distribution of 
these RRab, allowing us to trace its structure in a region 
that was still almost unexplored between b=-14°and b=- 
4° . We present the results of the fit of different models to 
the density profile of Sgr. The best fit to the core of Sgr is 
an exponential with a scale length of 4.1°along the major 
axis. When we look at the extension of Sgr we find a break 
(significant at the ~2er level) in the slope of the surface 
density along the main axis of Sgr. The nearly flat (or at 
least very slowly decreasing) profile in the outer region of 
Sgr shows that this dwarf galaxy is probably extending 
even further out our fields. 

Key words: Galaxies: dwarf - Galaxies: individual: Sagit- 
tarius dwarf - Local group - Stars: Variables: RR Lyr 



1. Introduction 

The Sagittarius dwarf galaxy is the closest known member 
of the Local Group orbiting around the Milky Way (~25 
kpc from the sun, ~16 kpc from the Galactic Centre), but 
as a consequence of its location behind the Galactic Cen- 
tre, it has been discovered only recently (Ibata, Gilmore, 
Irwin 1994, 1995). Since this discovery it turned out that 
Sgr presents typical features of a dwarf spheroidal: domi- 
nation of an o ld (> 10 Gyr) metal p oor st ellar population 
(Mat eo et al. |l995| ; Fahlm an et al. |l996| Marconi et al. 
1998 ; Bel lazzin i et al |l999| ) and absence of gas (Burton & 



Lockman 1999] ). Its highest surface density region is cen- 
tred on the Globular Cluster M54 (1=5.6°, b=-14.0°) and 
it is oriented roughly perpendicular to the Galactic plane 
so that its Northern extension (in Galactic coordinates) is 



completely hidden by the MW. 

The mapping of Sgr is difficult to achieve because of 
the combination of its low surface brightness (jiy > 
25.5 mag.arcsec -2 ), contamination by foreground Galactic 
stars and its large spatial extent (at least 22° x 8°) (Ibata 
et al. 1997, hereafter IWGIS). Evidence for the presence 
of Sgr has been established over 45°from b~ —3° (Alard 



1996] hereafter A96; Alcock et al 199 7, hereafter Alc97) 

down to b 48° (Mateo et al. |1998| , hereafter MOM), 

but it is difficult to assess whether these regions still cor- 
respond to the main body of Sgr or if we are merely en- 
countering tidal debris (as suggested by Johnston et al. 
1999). IWGIS proposed a map of the Southern part of 
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Sgr based on the spatial distribution of the bright main 
sequence stars in Sgr and covering an area of ~ 150 deg 2 
from b ~ —11° down to b ~ —26°. However, their method 
based on statistical decontamination fails at low Galactic 
latitudes (|b| <12°) where differential reddening and high 
density of foreground stars (only ~1 star in 1 000 is in Sgr 
in these regions) preclude any reliable decontamination, 
leaving the structure of the Northern extension of Sgr al- 
most unknown. To this point, the detection of RR Lyrae 
constitutes an essential tool to trace the structure of Sgr in 
these regions as they can be clearly separated from the RR 
Lyrae of the MW. This method has already proven suc- 
cessful and ~ 350 RRab were detected between b=-10°and 
b=-4°(A96; Alc97). However, a connection between these 
stars and the centre of Sgr was necessary in order to offer 
a clear vision of this important region strongly interacting 
with the MW. 

In this paper we report the detection of ^1500 RRab 
members of Sgr and located in its Northern extension. We 
present a surface density map of Sgr covering ~50 deg 2 
between b=-14°and b=-4°, based on the spatial distribu- 
tion of these variables. 

The paper is organized as follows : in section 2 we present 
our data (observations and reduction). Section 3 is de- 
voted to the description of the selection process of RR 
Lyrae stars as well as a study of its completeness. We 
then describe the structure of Sgr (section 4). Finally we 
summarize our results and conclude in section 5. 
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Table 1. Table 1. Observations 



Field 


DUO 


SAG 


Season 


1994 


1996 


Number of plates 


82 


69 


Field centre (l,b) 


(3.1°, -7.1°) 


(6.6°,- 


Emulsion 


III a J, III a F 


4415 


Filter 


GG385, RG630 BG12 


Limiting Magnitude 


Bj ~20.5 


V~20.0 



in.; 



2. Data 

2.1. Observations 

The data discussed in this paper consist of two sets of pho- 
tographic plates and films taken with the ESO lm Schmidt 
telescope at La Silla Observatory (see Table 1), each of 
them covering an area of ~ 25 deg 2 on the sky. 
The first set of plates was part of the DUO project aimed 
at detecting microlensing events towards the Galactic 
Bulge (Alard & Guibert |1997| , hereafter AG97). This field, 
centred on Galactic coordinates (1=3.1°, b=-7.1°), has al- 
ready been processed and presented in A96. The second 
set is new and includes 69 films centred on a field shifted 
towards the centre of density of the Sagittarius dwarf 
galaxy and slightly overlapping with the former (1=6.6°, 
b=-10.8°). Throughout the remainder of this paper, we 
will call the first field DUO field while the new field will 
be referred to as the SAG field. 



2.2. Data reduction 

The plates were scanned at CAI/Paris Observatory with 
the high speed microdensitometer MAMAQ (Machine Au- 
tomatique a Mesurer pour 1' Astronomie) , yielding images 
with a pixel size of 10//m (~ 0.6".). 

The photometric reduction has been performed with the 
software Extractor written by Alard. The process is as fol- 
lows: first a reference catalogue is extracted from a plate 
of good quality (seeing< 1"). For all the other plates, a 
new extraction is performed (implying a new detection of 
each object) and the new catalogue associated to the ref- 
erence catalogue. The light curves were built in this way 
plate by plate and stored in a database. For more details 
on the photometric reduction process see AG97. The fi- 
nal sample contains light curves for ~ 14. 10 6 stars in the 
DUO field and - 6.10 6 stars in the SAG field. 



2.3. Photometry 



1 MAMA ( jhttp: / /dsmama.obspm.fr| ) is operated by INSU 
(Institut National des Sciences de l'Univers) and Observatoire 
de Paris. 
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Fig. 1. Approximate bandpasses used in SAG (solid line) 
and DUO (dotted line). 



2.3.1. Calibration 

The DUO field has been calibrated with a CCD sequence 
taken at the ESO/Danish 1.5 m telescope at La Silla. The 
photometric system for this field is Bj = B — 0.28(B — V) 
(Blair & Gilmore 1982| ) . The Emulsion/Filter combina- 
tion was different for the SAG field and consisted of a 
Kodak Tech-Pan 4415 emulsion together with a BG12 
Filter. The Tech-Pan 4415 emulsion is an extremely fine- 
grained, high resolution film with a sensitivity extending 
to 0.69 /jm. For more infor mation s about the 4415 emul- 
sion, see Phillipps & Parker ( 1993 ) and references therein. 
We were not aware of any photometric relation published 
for the band used in SAG . The resulting bandpasses for 
both fields are shown on Fig. |[ The calibration in SAG 
has been performed with a sequence provided by Sara- 
jedini & Layden ( 1995| ), located 12" north of the glob- 
ular cluster M54 and consisting of 1638 stars calibrated 
in V and I bands. A polynomial fit to the instrumental 
magnitude of these stars yielded the photometric system 
B; = V + 1.47 (V — I), where B; stands for the magnitude 
in the color band used in SAG . The scatter about this 
relation is 0.17 mag. 



2.3.2. Correction for extinction 

The reddening has been estimated separately for each 
field. For the DUO field we used the well known prop- 
erty that the color of RR Lyrae stars at minimum magni- 
tude is approximately constant and depends only slightly 
on the period and the metallicity (Sturch, 1966] ). A red- 
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dening map has been estimated for this field by com- 
puting the mean colors (Bj-R) of RRab in small re- 
gions of 10'xlO'. The corrected magnitude is Bj = 
Bj - 2.84 E(B - R) (Wesselink |l987| ). For the SAG field, 
where no color information was available, we used the 



extinction map of Schlcgel et al. (1998, hereafter SFD) 
which provides reddening estimates with a precision of 
16% for |b| >10°. However, the SAG field extends to 
b^ -8°where, according to SFD, the reddening map 
might become inaccurate. Fr om th e relation E(V — R) = 
0.74 E(B - V) (Cardelli et al. |l989| ; hereafter CCM) we de- 
rive E(Bj - R) /E(B - V) = 1.46. This ratio in the overlap 
between DUO and SAG yields 1.32 ± 0.24, in reasonable 
agreement with the theoretical expectation, showing that 
even at the western edge of the SAG field the SFD map 
provides a satisfactory estimation for the extinction. As- 
suming E(V-I)=1.55 E(B-V) from CCM and a normal ex- 
tinction law A v =3.10 E(B-V) we obtain B io = B t - 5.38 
E(B-V) for the de-reddened magnitude in the SAG field. 




3. Detection of RRab 

3.1. The selection process 

We will describe here the selection process of RRab stars. 
For the sake of homogeneity, we reprocessed the stars of 
the DUO field, using the same selection criteria as for the 
SAG field. The search for RRab in DUO has been per- 
formed through the Bj band. 

A first selection was performed by calculating the x 2 
about the mean magnitude (Xmean) f° r each light curve. 
Stars with Xmean > ^ were then searched for periodic- 
ity. This cut should select all variables with an ampli- 
tude > 0.3 mag. A first estimate of the period was done 



with the string minimization method of Renson (1978). A 
more accurate period was then searched in a small win- 
dow spanning 0.1 day around the first estimate, using 
a multi -harm onic periodogram method (Schwarzenberg- 
Czerny 1996 ). The next step was to fit a Fourier series 
(with up to five harmonics) to the folded light curve: 

n<5 



Bi = A + A n cos(nujt + 0„) 



The x%t about the fitted light curve was then calculated 
and all the stars for which XraUo = Jxmean/Xju > 2 
have been selected as variable stars. At this step of the 
process the sample contained ~ 7 000 variables. 
The selection for RR Lyrae stars has been performed 
through the Fourier coefficients: for each variable we 
calculated the ratio of the amplitude of the first harmonic 
relative to the amplitude of the fundamental harmonic 
i?2i = A 2 /Ai, and the phase difference 02i = 02 — 20i. 
Fig.|| shows a plot of i? 2 i versus 02 1 for all stars satisfying 
Xmean > 8 an d Xratio > 2. Several clumps lie in this 
figure. The most obvious one is located at R21 ~ 0.45 



Fig. 2. R21 versus 02 1 for all variables with a large ampli- 
tude (xliean > 8 ) and a wel1 ntted n S nt CUrve (Xratio > 2). 

Also shown is the ellipse for selection of RRab. 



, (f>2i ~ 0.7. This clump corresponds to RR Lyrae stars 
of Bailey type ab (hereafter RRab). For lower values 
of i?2i {i.e.. for more symmetric light curves) we can 
distinguish two other clumps: one centred on (i?2i ~0.2, 
02i ~3.2) and a shallower one at (i?2i ~0.15, 02 1 ~ 1-75), 
corresponding respectively to contact binaries and RR 
Lyrae of Bailey type c (RRc). A faint strip across the 
plot at 021 ~ 3.1 is also visible and represents eclipsing 
binaries of Algol type. The selection of the RRab has 
been made with an ellipse centred on the clump (see Fig. 
|) and finally a cut on periods (0.40 rf > P > 0.85 d ) has 
been applied. The final sample contains ~ 3 000 RRab. 
The selected RRab may belong either to the MW or 
to the Sagittarius dwarf galaxy and we separated them 
through their distance modulus, assuming absolute 
magnitudes M Bj =0.79 (Wesselink |1987Q and My= 



=0.6 

(Mateo et al 1995). Furthermore, we take the mean color 
(V-I)o=0.46±0.06 after averaging over 27 RRab covering 
a wide range of metallicities from Table 1 of McNamara 
(1997). The apparent magnitude of each RRab has been 
estimated with the constant term of the Fourier series. 
Taking into account errors on the absolute magnitudes, 
apparent magnitudes, extinction and colors of RRab, 
the error on a single distance modulus is ^0.3 mag in 
both field, the main source of uncertainty coming from 
extinction. Fig|| shows the histogram of distance modulus 
for both fields before and after correction for extinction. 
The histograms were smoothed by estimating the mean 
magnitude every 0.1 mag in a 0.3 mag bin. Both his- 
tograms exhibit similar features: a broad bump centred 
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Fig. 4. RRab detected in the Sagit- 
tarius dwarf galaxy. The eastern box 
centred on (l,b)=(6.6°,-10.8°) repre- 
sents the SAG field and the west- 
ern box centred on (3.1°,-7.1°) is the 
DUO field. Each field covers an area 
of ~ 25 deg 2 . This map contains 
about 1500 RRab. The slight discon- 
tinuity in the density between DUO 
and SAG is due to the different com- 
pleteness levels of the plates (see 
text). Also shown are the RRab de- 
tected by the MACHO team (aster- 
isks). This map confirms that these 
RRab are the continuation of Sgr. 
Seven of their RRab are in common 
with ours. 



on (m-M)o ~14.5 (8 kpc) corresponding to RRab of the 
MW, and a sharp bump centred on (m-M)o ~16.9 (24 
kpc) representing RRab members of Sgr. According to 
current models of RRab densities in the Halo the Galactic 
contribution to the histograms for (m-M)o >16.3 should 
be no more than 5-10% (Wetterer & McGraw 1996). 
The 2D spatial distribution of all the RRab with a 
distance modulus greater than 16.3 is displayed on Fig.|j. 
This map includes ~ 1 500 RRab. The eastern and 
western box represents respectively the SAG field and the 
DUO field. The total area covered is about 50 deg 2 , and 
comprises the globular cluster M54 at (l=5.5°,b=-14.0°) 
which is associated to Sgr and located in its highest 
density region. The image of M54 is completely saturated 
until ~1.5 half mass radius on our plates, thus we do not 
expect this globular cluster to contribute significantly 
to our RRab sample. The spatial distribution of RRab 
reveals a density gradient in the SE-NW direction. 
We also show in Fig. |] the RRab discovered by the 
MACHO team (Alc97), confirming that these stars are 
the continuation of Sgr 



3.2. Completeness 

There are two steps where the completeness of the RRab 
sample might be affected: first, the detection of stars be- 
comes difficult towards the Galactic Centre because of 
the increasing stellar density, and some RRab blended by 
a neighbouring stars are missed. Second, we might miss 
some RRab during the selection process. 



3.2.1. Completeness of the extraction process 

To quantify the loss induced by the first effect we sim- 
ulated a set of 250000 artificial stars with the same ap- 
parent magnitude than the detected RRab. These stars 
were then injected in small regions of 10' x 10' uniformly 
spread over the fields and we tried to retrieve them with 
the same detection process as for the real stars. The lower 
panel of Fig.|| displays the fraction of stars re-detected 
as a function of Galactic latitude. Filled circles are stars 
injected onto the SAG field whereas crosses are stars sim- 
ulated in the DUO field. The dispersion reflects mainly 
the dependence on Galactic longitude. One can see that 
while the fraction of stars re-detected in SAG stays at a 
high level (above 90%) and varies slowly, this is not the 
case in DUO where this fraction drops abruptly to reach 
40% at b^-5°. The reason for the higher variation rate 
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Fig. 3. Upper panels (a): apparent magnitude histogram 
of the RRab stars. Lower panels (b): distance modulus 
histogram of the RRab stars. The dotted line indicates 
the distance modulus cut (16.3) adopted for membership 
of the Sagittarius dwarf galaxy. Note the reinforcement 
of the separation of the two bumps after correction for 
extinction. 



in DUO is that the stellar density gradient increases at a 
higher rate towards the Galactic Centre (the mean density 
gradient is ~15 stars. arcmin _2 .deg _1 in SAG and ~25 
stars. arcmin _2 .deg _1 in DUO ). Another feature visible 
on Fig.^b is that the loss induced by crowding is intrinsi- 
cally higher in DUO than in 571 G as can be seen in the 
range -10° <b<-8°(— 10% offset). This can be explained 
by the lower resolution of the III a j emulsion in DUO rela- 
tive to the finer grained 4415 emulsion in SAG (Parker & 
Malin 1999). Furthermore, the lower extinction in DUO 
^0.7) results in a higher number of stars de- 
S (SAG)~1.25 in the overlap), 



(A Bj /A Bi 

tected (N stars (DUO)/N 
increasing by this way the crowding. 



3.2.2. Completeness of the selection process 

Amplitudes: Our selection process might not be able to 
detect variable of low amplitude. To check the dependency 
of completeness on amplitude we simulated a set of 1 000 
RRab light curves with the same time sampling as the 
real ones, the Fourier coefficients have been taken from 
Simon & Teays (1982). The distributions in amplitude, 
magnitude and period (excluding integer fractions of a 
day) of the simulated light curves were chosen in a way 
to match the actual distributions of the detected RRab, 
and the phasing was uniformly distributed between and 



Fig. 5. Upper panel (a): detection efficiency for point 
sources as a function of amplitude. Each curve has been 
averaged over 50 000 simulated light curves. The full line 
corresponds to the SAG field and the dashed line cor- 
responds to the DUO field. The vertical dotted (resp. 
dash-dotted) line shows the amplitude cut adopted in 
DUO (resp. SAG ) for the construction of the surface 
density map. Lower panel (b): fraction of simulated stars 
rc-dctcctcd as a function of Galactic latitude. Filled cir- 
cles represent fields in SAG whereas crosses correspond to 
fields in DUO . 



2-7T. This set of simulated light curves was then injected 
in 100 regions of 10' x 10' (uniformly distributed over the 
fields) from which we took the errors to deteriorate the 
light curves. These light curves were then reduced in the 
same way as the real RRab. Fig.[i|a shows the completeness 
levels of our selection process as a function of amplitude 
for the two fields, averaged over 50 000 simulated light 
curves. The shapes of the completeness curves are nearly 
identical for both fields and the difference is not signif- 
icant. Fig.[|a shows that the detection rate stays above 
95% for amplitude >0.8 mag, and then drops abruptly 
down to ~ 20% at amplitude=0.5 mag. These results sig- 
nify that our selection process would detect almost all the 
RRab with an amplitude above 0.8 mag if these were point 
sources. 

However, the completeness levels will differ between 
SAG and DUO because the amplitudes of the RRab mea- 
sured in each filter are different, being more important on 
average for SAG than for DUO . This difference occurs 
because the color band of SAG peaks at shorter wave- 
length than the color band used for DUO whereas the 
amplitude of RRab decreases with increasing wavelength 
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(Smith 1995). A least square fit between the amplitudes 
of 30 RRab in common in the overlap yielded the relation 



A 



DUO 



0.98(±0.10)A 



SAG 



0.05(±0.11) 



(1) 



where Aouo and A-sag represent respectively the ampli- 
tudes measured in DUO and in SAG . In order to con- 
struct a consistent density map, we will consider in the 
remainder of this paper only those RRab satisfying am- 
plitude > 0.60 mag in SAG and amplitude > 0.54 mag 
in DUO, where 0.54 has been derived from the above re- 
lation (0). These cuts have been chosen both to ensure 
the largest sample as possible and to keep the complete- 
ness corrections at a manageable level. The corresponding 
corrections are 3.7% in SAG and 12.3% in DUO . 

Periods: Some RRab are missed because their periods 
are close to an integer fraction of a day, this causes points 
of the folded light curve to accumulate in a narrow phase 
range. The fitted Fourier series is then poorly constrained 
over a large fraction of the light curve and some of these 
stars might lie outside the ellipse of our selection process 
(see Fig.||). Monte-carlo simulations shows that we miss 
about ~ 30% of the RRab within the range 0.49 d to 
0.51 d for both fields, corresponding to a total loss of ~ 3%. 



18.54 18.53 

a }20oo (h° urs ) 



Fig. 6. Overlap region between DUO and SAG . The 
RRab detected in DUO (amplitude(Bj)>0.54) and SAG 
(amplitude(B J )>0.60) are shown respectively as crosses 
and open circles. The dashed lines represent the plate lim- 
its. Note that the scales are different in a and S. 



3.3. Homogeneity between the fields 

An important point to inspect for checking the consistency 
between the two fields is the overlap (Figra). Open circles 
correspond to RRab detected in the SAG field with an 
amplitude(Bi) > 0.60 mag. while crosses represent RRab 
detected in the DUO field and having an amplitude(Bj) 
> 0.54 mag. The dashed lines indicate the limit of each 
plate, their inclination is due to a slight tilt between the 
two plates. Note that this overlap applies to the reference 
frames and is not necessarily constant from plate to plate. 
Fig.|6| reveals that most of the RRab are detected indepen- 
dently in the SAG field and in the DUO field. However, 
some RRab are not detected twice and it is important to 
understand the reasons why these stars are missed by one 
of the fields: 

1. RRab not detected in the SAG field (crosses) 

— Eight RRab are located very close to the edge of the 
SAG plate. Such stars usually have fewer points in 
their light curves because the centre of the plate is 
not exactly the same at each exposure. For example 
the three westernmost stars in the SAG field have 
respectively 47, 55 and 47 points in their light curve 
(instead of 69 for most of the other light curves) . 

— Two RRab did not pass through the selection pro- 
cess (one because of its estimated Bj amplitude and 
the other one because of its Xratio)- 

2. RRab not detected in the DUO field (open circles) 



— One RRab has not been detected probably because 
of its low amplitude (0.61 mag in SAG ). 

— Two RRab have a period of nearly ~ 0.50 d and 
were detected in SAG only by chance. 

— Three RRab were blended by a nearby star. As 
stated above, DUO is more sensitive to crowding 
than SAG and the relative loss of three RRab in 
the overlap is fully consistent with the ~10% offset 
observed in FigJjpD in the range -10° <b<-8°. 

Most of the missed RRab will therefore have no statistical 
incidence and should not bias the density map. The only 
concern is for the greater sensitivity of the DUO field to 
crowding. However this effect should be lowered by the 
crowding correction. Turning now to the western edge of 
DUO we re-detect 7 RRab out of the 8 detected by the 
MACHO team in our field (disregarding two RRab located 
close to the edge). This is a satisfactory result. 

4. Structure of the Sagittarius dwarf galaxy 

4-1- Surface density of Sgr 

A surface density map is constructed from RRab with the 
amplitude cuts stated above. The spatial distribution of 
these RRab has been convolved with a Gaussian on a grid 
with a step of 0.1° and a variable filter size adapted to 
the local surface density a oc p^ 1 ^ 2 , constrained between 
0.2°and 0.5°. This map was then corrected for the different 
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a J3000 

Fig. 7. Smoothed map of the Sagittarius dwarf galaxy. This map is based upon the spatial distribution of RRab 
with distance modulus greater than 16.3. Only those RRab with (amplitude in B.;)>0.60 in SAG and (amplitude in 
Bj)>0.54 in DUO have been used. Completeness corrections have also been applied (see text). Contours are labelled 
as number of RRab per square degree. The dotted line (labelled 5) is not equidistant from the other contours. 



completeness in amplitude and crowding (see section 3.2). 
The resulting map is shown on Fig |7] where the elongated 
shape of Sgr is clearly visible. This is the first map of Sgr 
in these regions, showing that Sgr extends far beyond the 
outer limit of the map previously published by IWGIS. 
One of the most striking features of this map is the slow 
decrease (if any ?) of the density along the main axis of 
Sgr for |b| <9°. 

The main source of uncertainty in the surface density is 
the Poissonian noise in the star counts, which is variable 
over the field and tends to increase towards lower |b|. To 
estimate this noise we simulated 1 000 maps by injecting 
1 400 stars (corresponding to the number of RRab actually 
used to construct the final map) onto the field with a 
probability density matching the surface density of the 
real map. These spatial distributions were then processed 
exactly in the same manner as the real one and a la "noise 
map" has been deduced. This map is shown on Fig.^J 
where the contours are labelled in percent. The typical 



(relative) uncertainty is constrained between 10% and 15% 
over the main part of the field, but increases up to ^40% 
towards the edges where the number of RRab drops. 

4-2. Surface density profile of the main axis 

The position angle of Sgr has been determined by fitting 
an exponential to the surface density along various direc- 
tions. The highest scale length was reached for an angle 
of 108.4°, which we choosed as the direction of the major 
axis. Fig.^J displays the density profile of Sgr along that 
axis. This figure is based on a map smoothed on a constant 
scale of 0.5°. The thick line is the density after correction 
for completeness whereas the dotted line is the density be- 
fore that correction. The shaded region represents the la 
uncertainty issued from simulated maps. A discontinuity 
in the slope is clearly visible at ^6°from the centre. After 
this point the surface density seems to be almost constant. 
It is however disconcerting that this discontinuity occurs 
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Fig. 8. Uncertainty of the density map of Sgr based on 
lcr Poissonian error in the star counts. The contours are 
labelled in percentage. The dotted line (labelled 15) is not 
equidistant with the other contours. 




Distance from the center (degree) 

Fig. 9. Cross-section along the main axis of Sgr based on a 
density map smoothed with a constant filter size of 0.5°. 
The thick line represents the density after completeness 
corrections and the dotted line is the uncorrected density. 
The shaded region corresponds to the la uncertainty is- 
sued from simulated maps. Note the discontinuity in the 
density gradient at ~6°, also perceptible in the uncor- 
rected density. 



near the limit between DUO and SAG and we may won- 
der if this is not an experimental effect. We have shown in 
section 3 that the crowding correction in SAG and DUO 
were consistent with the completeness of each field in the 
overlap. Furthermore the break is also perceptible in the 
uncorrected density so it cannot be an effect of the crowd- 
ing correction. Another possibility is that our amplitude 
cut in DUO is to low to be consistent with SAG . This is 
difficult to check and we can only rely on those 30 RRab in 
common in the overlap, from which we derived the relation 



in amplitude between DUO and SAG (Eq. |l|). However, 
if we make the assumption that the RRab population is 
homogeneous over the field, it is possible to search for 
the relation Aouo = a, Asag + b for which the amplitude 
distributions are the most similar (through Kolmogorov- 
Smirnov test). The resulting coefficients were a=0.95 and 
b=-0.08 and are within the error bars stated in Eq. [j] 
The corresponding cuts are Asag = 0.6 <-> Aduo = 0.49. 
These cuts would have reinforced the discontinuity, show- 
ing that our adopted amplitude cuts are not responsible 
for the break observed in Fig. ^. We conclude that the dis- 
continuity in the slope of the density profile is probably 
real and not a consequence of the change of field. 
It is also possible to derive an upper limit for the extension 
of Sgr along the line of sight: the distance modulus his- 
togram can be roughly fitted by a Gaussian with a width 
of 0.2 mag, corresponding to a depth of ~4.5 kpc for an 
assumed distance of 24 kpc. 



Fig. 10. Location of the axes on which various models 
have been fitted. Each point is represented by the number 
of RRab (corrected for completeness) within a box of size 
0.5° x 1.0°. The vertical dashed line at a ^18.55 indicates 
the limit of the Maj6 axis (see text). 



4.3. Model fitting 

We define the following analytical functions to fit to the 
density profile: 



Pk = k 



1 



1 



[1 + (r/r c ) 2 ]V2 [1 + (r t /r c )2]V2 



Pe = p e e 



PG = Pg e- r2 ^ 



PL = a + a x r 



(2) 

(3) 

(4) 
(5) 
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Fig. 11. Results of the fit to the surface density of Sgr. Left panels: density profile along the main axis. The solid line 
is the fit to the whole main axis (MajlO) while the dotted line corresponds to the fit to the inner part of the main 
axis (Maj6). The dotted line is prolonged beyond 6°to allow visual comparison between the data and the inner fit. 
Right panel: density profile along an axis rotated 50°relative to the main axis (Min ). The errors bars represent the 
Poissonian noise. The fitted model is indicated in each panel. 



Where r represents the distance from the Centre of 
Sgr, and all other parameters are variables to be fitted. 
Eq. [2] refers to the empirical King model (K) with a core 
radius r c and tidal radius r t (King 1962). Eq. || refers to 
an exponential model (E) with a radius r e . Eq. ^ refers 
to a Gaussian (G) with a width of a g . Finally Eq. || 
refers to a linear model (L) where the density profile is 
modeled by a straight line. These models have been fitted 
along three segments. Two of these segments are located 
on the main axis: one corresponding to the main axis 
over its entire length (~10°), referred to as MajlO; and 
another one corresponding to the portion of the main axis 
contained within SAG (<6°from the center), referred to 
as Maj6. The latter segment has been chosen in order to 
avoid fitting the stars past the break and also because 
it is more consistent since it is entirely contained within 
SAG . Finally the minor axis is not present within our 
field and instead we fitted an axis making a large angle 



relative to the main axis (50 °), referred to as Min . The 
fit of model L on Majl Ohas been performed by fitting the 
density before 6°and after 6°separately. In order to get 
uncorrelated points for the fit, we took the densities (after 
correction for completeness) of RRab inside boxes with a 
size of 0.5° x 1.0°located along each axis (see Fig JlCj) . The 
results of the fit are shown in Table 2 and in Fig.|ll|. 
The single function model that best fits MajlOis model 
E (xj it /N£>oF—3-08)- However, the fit is significantly 
improved if we consider the model L ((xf it /NDOF—^-77) 
which reproduces the break already observed in Fig. |[ 
A Fisher test shows that the probability for the ratio of 
the Xfit/^DOF of these two fits to be lower than the 
observed value by chance is only ~13%. Note that we 
were unable to fit any convergent two-component model 
to MajlO: this is due to the almost constant density of 
the external region which causes one of the component 
to increase as we move away from the centre in order to 
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Table 2. Table 2. Results of the model fitting to the sur- 
face density of Sgr. The three first columns give respec- 
tively the axis, model and parameter considered. Column 
4 gives the result of the fit. Column 5 gives the uncertainty 
on the parameter value as given by the covariance matrix 
issued from the fitting procedure. Column 6 gives the x 2 
of the fit normalized by the number of degree of freedom. 
All values given in kpc assume a distance to Sgr of 24 kpc. 



Axis Models Parameters Value 


a 


X/«/ n dof 


MajlO K 


k 


178.32 


68.28 


4.11 




r c (kpc) 


2.69 


0.56 






r t (kpc) 


oo 






E 


pe 


120.87 


5.96 


3.08 




r e (kpc) 


2.42 


0.15 




G 


Pg 


84.53 


3.75 


5.10 




<j s (kpc) 


2.13 


0.08 




L 


.(in) 

a 


117.51 


5.64 


1.77 




-(in) 
a i 


-15.40 


1.38 






-(out) 
<Aq 


32.13 


11.33 






-(out) 

a l 


-0.50 


1.36 





Maj6 


K 


k 

r c (kpc) 
r t (kpc) 


139.91 99.88 2.59 
1.77 0.65 

oo 




E 


Pe 

r e (kpc) 


138.62 14.09 2.05 
1.73 0.21 




G 


Pg 

a s (kpc) 


104.43 5.04 2.46 
1.50 0.08 




L 


ao 
ai 


117.51 5.64 2.08 
-15.40 1.38 


Min 


K 


k 

r c (kpc) 
r t (kpc) 


344.84 108.04 2.00 
1.45 0.27 
2.98 0.30 




E 


Pe 

r e (kpc) 


160.12 11.79 5.03 
0.79 0.05 




G 


Pg 

a 9 (kpc) 


107.04 6.01 1.75 
0.84 0.03 




L 


ao 
ai 


99.30 4.69 3.16 
-18.97 1.00 



compensate the decrease of the other component. 
Concerning the core of Sgr (Maj6), the density 
profile is equally well fitted by model E and L 
(Xfit/NooF ~2.1). The scale length derived from 
model E is ~4.1°±0.5°(1.7±0.2 kpc). This value is 
slightly lower to the one derived by MOM who find an 



inner scale length of 4.7°in the Southern part of Sgr. 
Model K and G also give an acceptable fit to the core of 
Sgr but they fail to reproduce the high density in the first 
bin. Furthermore, the uncertainties on the parameters of 
the empirical King model are quite large and the infinite 
tidal radius is rather unrealistic. 

Finally, the best fit on Min is achieved by model G 
(Xf it /NDOF=l-75), but again it fails to reproduce the 
high density of the first bin. The only model that 
reproduces the high central density is model E but the 
Xfu/^DOF of this model is worsened by the poor fit on 
the three last bins. However these bins contain only very 
few points (between 1 and 5), introducing uncertainties 
induced by small-numbers statistics. 



5. Conclusion 

To summarize, we presented the detection of ~1 500 
RRab stars located in the Sgr dwarf galaxy. A surface 
density map based on the spatial distribution of these 
variables unveiled the structure of this dwarf galaxy in 
a region that was still almost unexplored so far between 
b=-14°and b=-4°. The core of Sgr is best fitted by an 
exponential with a scale length of 4.1°along the major 
axis. A cross section of this density map revealed a break 
in the slope occurring at ~6°from the highest density 
region of Sgr and an almost flat density past the break. 
Although the break coincided with the change of field we 
have shown that this is unlikely to be an experimental ef- 
fect since it is also perceptible in the uncorrected density, 
whereas the DUO field is intrinsically more sensitive to 
crowding than SAG (lower resolution, lower extinction). 
Also, as shown in Section 4.2, the amplitude cuts used 
in this study cannot be considered as responsible for the 
break. Finally, could this break be a consequence of an 
overestimation of the completeness correction in DUO 
relative to SAG ? Though not excluded, this would be 
in conflict with what is observed in the overlap where 
3 RRab blended by a neighbouring star were detected 
in SAG and missed in DUO , a result that is quite 
consistent with the corrections actually applied. We argue 
therefore that the break is real. The significance of the 
break relative to an exponential with a scale length of 
4.1°is ~2<7. MOM also observed a break in their density 
profile in the Southern extension of Sgr. However neither 
the location (20°from the centre) nor the density at the 
break location (Xy ^29.0 mag.arcsec -2 ) are consistent 
with our values (6°and ~26.7 mag.arcsec -2 ) implying 
that cither the main body of Sgr is not symmetric or 
these "post- break" stars are not directly related to it. 
Another striking feature revealed by the surface density 
profile is its flatness past the break. This feature relies 
on the accuracy of the completeness correction over the 
field, a correction that becomes quite important at low 
Galactic latitudes (up to 60%). Yet, the difficulty of 
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modeling point spread functions on photographic plates 
(due to non-linear response of the emulsion) and potential 
systematic errors caused by differential sensitivity over 
the plate makes the crowding correction rather uncertain. 
Therefore, although our completeness corrections are 
fairly consistent within the overlap, we cannot exclude 
that the flatness of the density profile in the outer regions 
is a consequence of an overcorrection. Wide-field high 
resolution imaging would be necessary in these extremely 
crowded regions (up to ~10 6 stars per square degree at 
our magnitude limit) to confirm or to rule out this issue. 
Nevertheless, even if we consider that our completeness 
corrections are overestimated by a factor of 2 (a quite 
conservative estimate), it remains that the density profile 
decreases slowly in the outer regions and Sgr may well 
be extending even further out towards (beyond ?) the 
Galactic plane. 



Johnston et al. ( 1999 ) recently modeled the Sgr stream as 
a superposition of a main body and tidal streams of stars 
stripped on previous peri-centric passages. This scenario 
has been worked out to explain both the break observed 
by MOM and the possible detection of stars in the outer 
region of Sgr with different radial velocities relative to 
those of the main body (Majewski et al. 1999). Similarly, 
spectroscopic observations on our RR Lyrae catalogue 
could allow to determine the nature of the stars in the 
outer region: if these stars are linked to the main body 
of Sgr, then they should share almost the same radial 
velocities as the main body (apart of a gradient along the 
main axis due to the rapidly varying Galactic potential). 
On the other hand, if the break we observe corresponds 
to a transition between the main body and an unbound 
tidal stream from a previous orbit, it is likely that the 
two objects will have different radial velocities. This 
new catalogue of RR Lyrae is an interesting opportunity 
to study further a region of Sgr that has been poorly 
investigated so far. 
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